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Abstract 

We address a decentralized wireless communication network with a fixed number u of frequency sub-bands 

to be shared among N transmitter-receiver pairs. It is assumed that the number of users N is a random variable 

with a given distribution and the channel gains are quasi-static Rayleigh fading. The transmitters are assumed to be 

^ unaware of the number of active users in the network as well as the channel gains and not capable of detecting the 

CO 

presence of other users in a given frequency sub-band. Moreover, the users are unaware of each other's codebooks 

in 

and hence, no multiuser detection is possible. We consider a randomized Frequency Hopping (FH) scheme in which 
each transmitter randomly hops over a subset of the u sub-bands from transmission to transmission. Developing a 
new upper bound on the differential entropy of a mixed Gaussian random vector and using entropy power inequality, 
we offer a series of lower bounds on the achievable rate of each user. Thereafter, we obtain lower bounds on the 

•l-H 

^ maximum transmission rate per user to ensure a specified outage probability at a given Signal-to-Noise Ratio (SNR) 

H 

level. We demonstrate that the so-called outage capacity can be considerably higher in the FH scheme than in the 
Frequency Division (FD) scenario for reasonable distributions on the number of active users. This guarantees a 
higher spectral efficiency in FH compared to FD. 
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I. Introduction 

A. Motivation and Related Works 

Optimal resource allocation is an imperative issue in wireless networks. Wide applications of wireless 
systems, in recent years, and the limited available resources in the network necessitate efficient usage 
of such resources. Multiuser interference is known to be the most important factor, which degrades the 
network performance when multiple users share the same spectrum. There has been a tremendous amount 
of research on designing an efficient and low complexity resource allocation scheme that maximizes the 
quality of service per user while controlling the detrimental effect of multi-user interference. Existing 
resource allocation schemes in the literature are classified as centralized, i.e., a central controller manages 
the resources, or decentralized, where resource allocation is performed locally at each node. 

The main goal in traditional wireless systems was to avoid the interference between users by transmitting 
over orthogonal channels. A well-known example of such systems is the Frequency Division (FD) system, 
in which different users transmit over disjoint frequency sub-bands. The assignment of frequency sub-bands 
is performed by a central controller. Despite its simplicity, [1] proves that in a wireless network in which 
the interference is treated as noise (no multi-user detection is performed), and if the crossover gains are 
sufficiently greater than the forward gains, FD is Pareto-optimal. Due to practical considerations, such FD 
systems have a fixed infrastructure, i.e., they rely on a fixed number of frequency sub-bands. If the number 
of users changes, the system is not guaranteed to offer the best possible spectral efficiency because, most 
of the time, the majority of the potential users may be inactive. In recent years, many centralized power 
and spectrum allocation schemes have been extensively studied in cellular and multihop wireless networks 
[2]-[9]. Clearly, centralized resource allocation schemes provide a significant improvement in the network 
throughput over decentralized (distributed) approaches. However, they require extensive knowledge of 
the network configuration. In particular, when the number of nodes is large, deploying such centralized 
schemes may not be practically feasible. 

Most of the decentralized schemes in the literature rely on either game-theoretic approaches or cognitive 
radios. Cognitive radios [10] have the ability to sense the unoccupied portion of the available spectrum. 
Fundamental limits of wireless networks with cognitive radios are studied in [1 1]— [15]. These smart radios 
require sophisticated techniques do detect the spectrum holes that add to the overall system complexity 
[16]. As such, it is essential to have a decentralized spectrum sharing strategy without using cognitive 
radios, which allows the users to coexist while utilizing the spectrum efficiently and fairly. 
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Due to its interference avoidance nature, hopping is the simplest spectrum sharing method to use in 
decentralized networks. Frequency Hopping (FH) is the most popular scenario in this category in which 
users randomly switch to different frequency sub-bands from transmission to transmission. As different 
users typically have no prior information about the codebooks of the other users, the most efficient method 
is avoiding interference by choosing unused channels. As mentioned earlier, searching the spectrum to 
find spectrum holes is not an easy task due to the dynamic spectrum usage. As such, FH is a realization 
of transmission without sensing, while avoiding the collisions as much as possible. Frequency Hopping is 
one of the standard signaling schemes [17] adopted in ad hoc networks. In short range scenarios, bluetooth 
systems [18]-[20] are the most popular examples of a wireless personal area network, or WPAN. Using 
FH over the unlicensed ISM band, a bluetooth system provides robust communication to unpredictable 
sources of interference. A modification of Frequency Hopping called Dynamic Frequency Hopping (DFH) 
selects the hopping pattern based on interference measurements in order to avoid dominant interferers. 
The performance of a DFH scheme when applied to a cellular system is assessed in [21]-[23]. 

Although there has been a tremendous amount of work on the performance evaluation of hopping-based 
decentralized networks, there are only a few information-theoretic results reflecting the fundamental limits 
of such networks. In a pioneering work, [24] offers a clean analysis of a decentralized network where all 
users follow a randomized hopping strategy to share the resources, while the number of active users is a 
random variable with a given distribution. In [24], the channel gains and the number of active users are 
assumed to be static and known to the corresponding transmitter-receiver nodes. This assumption makes 
the concept of achievable rate in the Shannon sense meaningful. However, in case the channel gains and 
the number of active users are not known to the transmitters, the concept of achievable rate may no longer 
be valid. A common setup for such an assumption is a network where channel gains are quasi-static fading 
and unknown to the transmitters, which is the framework for this paper. 

Rayleigh fading is an unavoidable phenomenon in wireless networks that can affect the performance of 
the system significantly. Traditionally, Rayleigh fading has been considered to be harmful due to reducing 
the transmission reliability in wireless networks. However, recently, researchers have been able to reduce 
this harmful effect by exploiting the so-called multiuser diversity [25], [26]. This can be considered a 
scheduling gain by allowing the users with favorable channels to be active. It is shown that multiuser 
diversity gain can be as large as log log N in broadcast and multiple-access channels [27]-[29] and logiV 
in single-hop ad hoc networks [30], [31], as N grows to infinity. However, achieving this scheduling gain 
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requires the fading channels to vary over time such that all possible realizations of the fading process are 
covered. In the case that channel gains are selected randomly at the start of the transmission and remain 
constant during the whole transmission period (quasi-static fading), the channels do not have ergodic 
behavior. In this case, a suitable performance measure is the e-outage capacity [32], denoted by R(e), 
which is defined as the maximum transmission rate per user, ensuring an outage probability below e, i.e., 

R(e) = sup{R : Pr{Outage} < e}. 

The reality of wireless channel is more complicated to be simply represented by Rayleigh fading model. 
A class of channel models considered in the literature is the one in which the signal power decays according 
to a distance-based attenuation law [33]-[40]. Moreover, the presence of obstacles adds some randomness 
(known as shadowing) to the received signal. It is well known that the effects of such random phenomena 
can significantly affect the throughput of a spectrum sharing network in both multi-hop [41]-[44] and 
single-hop scenarios [45] (Chapter 8), [46]-[50]. These features indeed increase the frequency reuse factor 
as they will attenuate the interference caused by a given transmitter on its neighboring receivers. In spite of 
the significance of the effects of distance-based attenuation and shadowing on the throughput of a spectrum 
sharing system, unfortunately, there is not a single commonly accepted model to represent these factors. 
It should be emphasized that the inclusion of signal attenuation due to distance and/or shadowing will 
indeed simplify the spectrum sharing as the multi-user interference will be attenuated and consequently 
its harmful effect will be reduced. On the other hand, these factors do not impact the performance of the 
orthogonal schemes in which the multi-user interference is altogether avoided. In spite of this fact, as the 
actual model used to represent distance-based attenuation and shadowing can have a profound impact on 
the system throughput (to the advantage of the randomized spectrum sharing schemes advocated in the 
current article), to avid any confusion, we have relied on a simple Rayleigh fading model which in some 
sense captures the minimum advantage offered by the proposed scheme vs. those based on orthogonal 
separation of users. 

Reference [51] studies a wireless network composed of a set of transmitter/receiver pairs in which a 
given link can be off or transmit with a constant power. [51] considers both the case of Rayleigh fading 
as well as a Rayleigh fading mixed with a proper distance-based attenuation and among other results 
provides a comparison between the scaling (with respect to the number of links) of the throughput in 
these two cases. 
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In [52], [53], the authors study a decentralized wireless ad hoc network where different transmitters 
are connected to different receivers through channels with a similar path loss exponent. Assuming the 
transmitters are scattered over the two dimensional plane according to a Poisson point process, a fixed 
bandwidth is partitioned into a certain number of sub-bands, such that the so-called transmission intensity 
in the network is maximized, while the probability of outage per user is below a certain threshold [52]. The 
transmission strategy is based on choosing one sub-band randomly per transmission, which is a special 
case of Frequency Hopping. In [53], a non-iterative and distributed power control scheme is introduced, for 
which the constant power and the channel inversion schemes are extreme cases. It is observed that none of 
these cases are ideal in general. In fact, it is shown that regulating the transmission power proportionately 
to the inverse square root of the forward channel strength minimizes the outage probability. 

Recently, Orthogonal Frequency Division Multiplexing (OFDM) has been considered as a promising 
technique in many wireless technologies. OFDM partitions a wide-band channel to a group of narrow-band 
orthogonal sub-channels. This motivates us to consider the underlying system to consist of u narrow-band 
orthogonal frequency sub-bands. 

B. Our Contribution 

In this paper, we consider a decentralized wireless communication network with a fixed number u of 
frequency sub-bands to be shared among N transmitter-receiver pairs. Any transmitter is connected to any 
receiver through a channel with quasi-static and non-frequency selective Rayleigh fading. In other words, 
the channel gains are picked randomly (based on Rayleigh distribution) at the start of the transmission and 
remain fixed for the whole transmission. It is assumed that the number of active users is a random variable 
with a given probability mass function. The channel gains and the number of active users are unknown 
to all transmitters, however, the receivers are assumed to be aware of their direct channel gains and the 
interference Probability Density Function (PDF). Moreover, users are unaware of each other's codebooks 
and hence, no multiuser detection is possible. A randomized Frequency Hopping scheme is proposed 
in which each transmitter randomly hops over v out of u sub-bands from transmission to transmission. 
Assuming i.i.d. Gaussian signals are transmitted over the chosen sub-bands, the distribution of the noise 
plus interference becomes mixed Gaussian, which makes calculation of the achievable rate complicated. 
The main contributions of the paper are: 

• Developing a new upper bound on the differential entropy of a class of mixed Gaussian random 
vectors and using entropy power inequality, we offer three lower bounds on the e-outage capacity 
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of each user denoted by -R F Hib( e )' -^FHib( e )' an( ^ ^FHib( e )- To evaluate the system performance 
analytically, we use -R F Hib( e )> which can be computed easily. However, computation of -R F Hib( e ) 

(2) 

and -RpHib( e ) involves integrations that cannot be carried out in a closed form. In the simulation 
results, we use the lower bounds -R F H,ib( e ) an( ^ -^FH,ib( e )' which are tighter than -R F H,ib( e )- 

• We perform asymptotic analysis for the outage capacity in terms of e and SNR. In the asymptotically 
small e regime, we observe that the maximum of outage capacity is obtained for either v — 1 or v — u. 
In the asymptotically small SNR regime, we demonstrate that for any value of v the system achieves 
the optimal performance. For asymptotically large values of SNR, it is shown that v opt = , 
where 

Kmax is the maximum possible number of concurrently active users in the network. 

• We compare the outage capacity of the underlying FH scenario with that of the FD scheme for 
various scenarios in terms of distributions on the number of active users, SNR and e. It is shown 
that FH outperforms FD in terms of outage capacity in many cases. We observe that in the low SNR 
regime, FH and FD offer the same performance. In the low e regime, FD is always better than FH, 
however, for many practical scenarios there exists a threshold, e t h, such that FH outperforms FD as 
far as e > e th . Also, we have shown that supremacy of FH over FD in the high SNR regime occurs 
quite often. 

The paper outline is as follows. The system model is given in section II. Section III describes analysis 
of the outage capacity. Section IV is devoted to derive lower bounds on the achievable rates of users. Also, 
in this section, we offer a new computable upper bound on the differential entropy of a mixed Gaussian 
random vector. In section V, based on the results in sections III and IV, we discuss how the users in the 
FH system fairly share the spectrum while maximizing the outage capacity. Derivation of various lower 
bounds on the outage capacity of users is part of the materials in this section. Comparison between the 
FH and FD scenarios is given in section VI through the analysis and simulation results. Finally, section 
VII concludes the paper. 

C. Notation 

Throughout the paper, we use the notation E{.} for the expectation operator. For a function f(X,Y) 
of two independent random variables X and Y, E x {f(X,Y)} denotes the expectation of f(X,Y) with 
respect to X while Y is treated as a parameter. Pr{£} denotes the probability of an event £, 1(8) the 
indicator function of an event 8 and px(-) the PDF of a random variable X. Also, l(X; Y) denotes the 
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mutual information between random variables X and Y and h(X) denotes the differential entropy of a 
continuous random variable X. 



II. System Model and Assumptions 

We consider a wireless network with N users^] operating on a bandwidth consisting of u sub-bands. It 
is assumed that the i th user exploits Uj out of the u sub-bands in each transmission and hops randomly to 
another set of Vi frequency sub-bands in the next transmission. This user transmits independent complex 
Gaussian signah 2 of variance — over each of the chosen sub-bands in which P denotes the total average 
power of each transmitter. Each receiver is assumed to know the hopping pattern of its affiliated transmitter. 
It is assumed that the users are not aware of each other's codebooks and hence, no interference cancellation 
is performed at the receiver sides. The quasi-static and non frequency- selective fading coefficient of the 
channel from the i th transmitter to the j th receiver is shown by hij. All the channel coefficients in the 
network are assumed to be complex zero-mean Gaussian random variables of unit variance corresponding 
to Rayleigh fading. However, due to the absence of any feedback link, the transmitters do not have 
information about any of the channel gains. By the same token, the transmitters are not aware of the 
number of active users in the network. 

As all users hop over different portions of the spectrum from transmission to transmission, no receiver is 
assumed to be capable of tracking the instantaneous interference. This assumption makes the interference 
plus noise PDF at the receiver side of each user be mixed Gaussian. In fact, depending on different choices 
the other users make to select the frequency sub-bands and values of the crossover gains, the interference 
on each frequency sub-band at the receiver side of any user has up to 2 N ~ 1 power levels^} The vector 
consisting of the received signals on the frequency sub-bands at the i th receiver in a typical transmission 
slot is 

Yi = hi : {Xi + Z{, (1) 

where Xi is the u x 1 transmitted vector and Zi is the noise plus interference vector on the receiver side 
of the % th user. Due to the fact that each transmitter hops randomly from transmission to transmission, 



Each user consists of a transmitter-receiver pair. 



2 Since in this work we deal with fading channels, the transmitted signals are assumed to be complex for simplicity of analysis. 
3 The knowledge of the channel gains and the number of active users at the receiver side is realized by identifying the interference PDE 
4 It is notable that the interference plus noise PDF has 2 iv ~ 1 power levels, almost surely, as the channel gains are considered to be 
continuous random variables. 
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one may write (■) as 

Px^) = Yl 7^\9u(x,C), (2) 
ce<t W 

which corresponds to the mixed Gaussian distribution. In the above equation, g u (x,C) denotes the PDF 
of a u x 1 complex zero-mean jointly Gaussian vector of covariance matrix C and the set <£ includes 
all u x u diagonal matrices in which Vi out of the u diagonal elements are ^ and the rest are zeros. 
Denoting the noise plus interference on the j th sub-band at the receiver side of the i th user by Z i: j (the 
j th component of Z,j), it is clear that p Zi ■ (•) is not dependent on j. This is by the fact that crossover gains 
are not sensitive to frequency and there is no particular interest to a specific frequency sub-band by any 
user. We assume there are + 1 (L, < 2 Ar ~ 1 — 1) possible nonzero power levels for Zij, say {of j}£i . 
Denoting the occurrence probability of a 2 t by a i:l , Pz i:i (-) is given by 

P*A Z ) = J2 ^ exp ( "S" ) ' (3) 

1=0 V ^ J 

where a 2 = af fi < of x < af 2 < ■■■ < of^. (o -2 is the ambient noise power). We notice that for each / > 0, 
there exists a > such that of ; = a 2 + q^P where = q i0 < Q,i < Q,2 < ••• < Q In fact, one may 
write = J2k=^ e k,jhk,iX k j + z/jj where X k j is the signal of the k th user sent on the j th sub-band, 
€k,j is a Bernoulli random variable showing if the k th user has utilized the j th sub-band and i/jj is the 
ambient noise, which is a zero-mean complex Gaussian random variable with variance a 2 . The ratio ^ 
is taken as a measure of SNR as is denoted by 7 throughout the paper. 

Since the transmitters are not aware of the channel gains and the number of active users in the network, 
the Shannon capacity is not meaningful in this setup. In this case, a suitable performance measure is the 
e-outage capacity, denoted by R(e), which is defined as the maximum transmission rate per user ensuring 
an outage probability below e, i.e., 

R(e) = sup{R : Pr{Outage} < e}. (4) 

III. Analysis of the Outage Capacity 

Let hi contain the channel coefficients concerning the i th user, i.e., hi = (hn ■■■ ft,^^ . In this 
case, we denote the achievable rate of the i th user by Mi(hi). The outage event for this user is 

Oi(R)±{hi-.Mi(hi)<R}, (5) 
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where R is the actual transmission rate. Hence, 

R(e) =sup{# : Pr{Oi(R)} < e}. (6) 

We emphasize that the randomness of the number of active users is involved in the outage event, as N 
represents the size of hi. Moreover, due to symmetry, the e-outage capacity is the same for all users. 

Having hi fixed, it can be observed that the communication channel of the i th user is a channel with 
state Si, the hopping pattern of the i th user, which is independently changing over different transmissions 
and known to both the transmitter and receiver. Hence, 

= KXfMlSi) = Pr (^ = s l )l(X l ;Y l \S l = 8i ), (7) 

where I(Xj; Yi\Si = si) is the mutual information between Xj and % for the specific sub-band selection 
dictated by Si = Sj. The set (5, ; denotes all possible selections of v i: out of the u sub-bands. As pg.(-) ls 
a symmetric density function, meaning all its components have the same PDF given in ([3]), we deduce 
that l(Xi]Yi\Si = Si) is independent of Sj. Therefore, to calculate &i(hi), we may assume any specific 
sub-band selection for the i th user in (3,, say the first Vi sub-bands. Denoting this specific state by s*, we 
get 

M l (h i ) = l(Xi;Y i \S i = s*). (8) 
In this case, we denote Yj, and X; by Yi(s*) and Xi(s*) respectively. Obviously, 

= KMs*); Y t (s*)) = h(y,«)) - h(Z0- (9) 

As Yi(s*) and Z { are complex mixed Gaussian vectors, there is no closed expression for the differential 
entropy of these vectors. As such, we provide a lower bound &ijb(hi) on &i(hi) in the following section. 
Subsequently, using &i y vo{hi), we derive a lower bound on the outage capacity of the i th user as 

R lh (e) = supji? : Pr{hi : ^,i b (^) < R} < e}, (10) 

and show that this lower bound is higher than the actual outage capacity in the FD scheme in many 
scenarios. 
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IV. Lower Bounds on Mi(hi) 

The aim of this section is to find a lower bound on £%i(hi). The idea behind deriving this lower bound 
is to invoke entropy power inequality (EPI). As we will see, this initial lower bound is not in a closed 
form as it depends on the differential entropy of a mixed Gaussian random variable. In appendix A, we 
obtain an appropriate upper bound on such an entropy, which leads us to the final lower bound on Mi(hj). 

Theorem 1 There exists a lower bound on £%i{h^) given by 



i(hi) > vdogi Li , 1 +1), (ID 



2-^2 y -|fe M | 2 7 

where Jg? = - £fio a u log a hl and % = ^ Efii Hi log (l + 

Proof: Let us define X[ to be the Vi x 1 signal vectoij^] of the i th transmitter that is sent through the 
first v i frequency sub-bands. Let Y? = h^X^ + Z[, where Z[ is the noise plus interference vector at the 
receiver side of the i th user on the first Vi frequency sub-bands. According to the classical EP^} 

2 v i 1 > 2 v i 1 + 2 v i 1 . (12) 

— h(Z') 

Dividing both sides by 2 fl < * , 

h(f/) - h(Zl) > Vi log ^M hi ^'^) + ij . (13) 
On the other hand, since Y( is a subvector of Yi{s*), 

= i(x t ( s ;y, Yi(s*)) > i(x; ; y>) = h(y/) - h(zQ. (14) 



Based on (13) and (14), we get the following lower bound on 3$i{hj), 

MSi) > ^ log(2^( h(fti - i ^ ) - h( ^ ) ) + 1). (15) 

Clearly, h(/ij jX-) = log (ne ^ hl '^ P j . As Z[ has a mixed Gaussian distribution, there is no closed formula 
for h(Z l / ). To circumvent this difficulty, we have to find an appropriate upper bound on h(Z^). 
We start with the following Lemma: 



5 X[ consists of the first Vi elements of Xj(sf). 

6 As we deal with complex random vectors, the format of EPI is different from its real counterpart. 
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Lemma 1 Let be a t x 1 complex mixed Gaussian random vector with different covariance matrices 
{Qih}i=i and corresponding probabilities {pi\\ =1 where q\ < ■ ■ ■ < q\ and J2t=iPi = 1- Then, 

L 

h(©) < tj^p/ log (Tree?) +^-^ (16) 
1=1 

where 

L 

^ = -J2p i1 °zp i (17) 

1=1 

and 

& = € x>i°g f 1 + 4^^) . (is) 

Proof: See appendix A. ■ 
Using the chain rule for differential entropy, 

Vi 

h(^)<5^h(Zy). (19) 

J'=l 

Applying the special case of Lemma [T] corresponding to t — 1 for a scalar complex mixed Gaussian 
random variable given in ([3]), 

Li 

HZ itj ) < Hi log (20) 



where 



-^2a i; iloga i; i (21) 



and 



a 2 U ( „2 s^l-1 



E ay log (1 + %^)^ 
1 Va^log f 1+ (^ + l)E^o«v, 



(22) 



Substituting (20) in (19) gives an upper bound on h(Z 4 ') as 

Li 



h(4) < E a ^ lo § (^J*) + u * - ^) • (23) 



Z=0 
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By d23) and (fBL 



i(hi) > Vilog \2 Vi \ V "* J * ,; V+l 



2- J *2*|fe M | 2 7 

.^nf4(w+i) a< 



^log( / 1 +1). (24) 



Let us define 



^ nf=i( c M7+i) 



- Vi lo s I : tt*« ^ J, iU., + 1 1 • (25) 



The above proof reveals the following observations: 

Observation 1- As Li = 2 N ~ l — 1 with a probability of 1, it can be immediately verified that J£? does 
not depend on the crossover gains. However, ^ is implicitly a function of all crossover gains as the partial 
sums J2m=i a i,m for 2 < I < L, depend on the ordering of the crossover gains. This will be investigated 
more in Lemma |2] 

Observation 2- Since Yltli( c i,n + ^) ai ' 1 — ll£i( c i,£ l 7 + l) a<,! = ( c i,L,7 + l) ( ' 1 ~ ai ' ' ) > one obtains a looser 
version of &fi h (h>i) given by 

We note that &f\ h {hi) still has the same asymptotic expression as that of &^\ h (hi) in the high SNR regime. 
As we will see later, the computaion complexity of the lower bound on the outage capacity inspired by 
&-%(hi) is much lower than that of M^ h {hi). 

Let us consider a "fair" FH system in which vi = v for some 1 < v < u and for any 1 < i < N. As 
explained before, Sfj depends on the ordering of {cy}^! , which requires analyzing the order statistics of 
the channel gains. To avoid this, the following Lemma introduces a lower bound on % that only depends 
on c i: Li i.e., the largest interference crossover gain. 

Lemma 2 In a fair FH system, 

*=-(W-l)(llogH + (l-^)log(l-^)) (27) 



and 



, > % fc A E * ( lQ S (1 + (1 ~ <w) } ~ (N - l)g log | 

4 ' lb Ci )Li 7 + 1 
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where B is a Binomial random variable with parameters [N — 1, 

Proof: Let us define p = ^. As Li = 2 N ~ 1 — 1 with probability one and each user selects a certain 
frequency sub-band with probability p. The collection {a^}^! consists of the numbers p>(l — p) N ~ 1 ~ j 
repeated ( iV ~ 1 ) times for < j < N — 1. Hence, 

N-l 



= -(e^'^J 1 )^ 1 -^" 1 "'') 10 " 

- ( f^v - 1 - j) - x y (i - p)^- 1 -^ io g (i - p) 

= -(TV - l)p logp - (TV - 1 - (JV - l)p) log(l - p) 

= -(7V-l)(plogp + (l-p)log(l-p)). (29) 

As for computation of Yl l m =G a i,m is not an easy task. In fact, it depends on the ordering of the 
crossover gains. For example, if TV = 4, 

3 r (1 _ p) 3 + 2p(1 _ p) 2 +p2(1 _ p) if 1^)2 <|/ l3)1 |2< l^p + l^p < |^ 4ji |2 

m =0 { (1 -p) 3 + 3p(l -p) 2 if |/l 2>1 | 2 < l^lP < Kl| 2 < l^2,l| 2 + |/>3,l| 2 

(30) 

To avoid this difficulty in describing 5fj, we derive a lower bound on this quantity, which is not sensitive 
to the ordering of crossover gains. Taking each a iy i, there exists aO<s<A r — 1 such that = 



J2 S _ \h k i\ 2 

p s (l — p)-^ -1-5 . This implies that a i){ corresponds to the interference plus noise power level 3=1 j ' 1 P + 

y^s |2 

a for some 1 < h < ■ ■ ■ < k s < N where kj ^ i for 1 < j < s. Since f^LL*^ p + a 2 > 
Ete ^ {1 - 2 ;-- s}|fefct ' t|2 p + a 2 for any set A C {1, 2, • • • , s}, and Ete ^ {1 ' 2 ;-- s} '^"'V + a 2 is itself a power 
level in the PDF of the noise plus interference on each frequency sub-band, we conclude that its associated 
probability p^(l — p)^- 1 -!- 4 ! is an element in the sequence (a i)0 , a i: \, ■ ■ ■ ,a i: i-i). Therefore, we come 
up with the following lower bound, 

l-l s-1 



m=0 _4C{l,2,-,s} s'=0 ^ ' 



(3D 
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Using (31) in (22) yields 

AT-l 



p.- > - rt-- * (> + ( ^ +i)E y y,(i - ri "" 

t + 1 V s / 



AT-l 



1 /AT _ 1\ 

— T E W - P)"~ 1_ ' lo S ( X + (! " P")^7) 



E B {log (1 + (1 - p B )c hLl l) }-{N- l)plogp 



(32) 



Cj.iiT + 1 

where B is a Binomial random variable of parameters (N — l,p). ■ 
From now on, we replace Sfj with Sfjjb in all expressions offered for the lower bounds on 3?i(hi). In 
a fair FH system, we denote a it0 , and % )Vo by a(v,N), Jt?(v,N) and ^ iilb (v,iV) respectively to 
emphasize their dependence on the parameters v, N. 

As a special case of the fair system, let us assume Vi = u for all i, i.e., all users spread their power 
on the whole spectrum. This scheme is called full-band spreading (FBS). In this case, it can be observed 
that a^i = for / < L { - 1 and a iM = 1. This yields a(u,N) = J>tf{u,N) = % b (u, iV) = 0. In fact, 
Sif^JJki) is tight for v — u, i.e., &f\ h (hi) is exactly the achievable rate of the i th user while all users 
transmit over the whole spectrum. We denote this rate by &i t FBs{hi), which is given by 

&i,FBs(hi) = ulog [ ifl'.7 7 -v + 1 ) • (33) 

Observation 3- A straightforward method to develop a lower bound on the achievable rate of an additive 
non-Gaussian noise channel is to replace the noise with a Gaussian noise of the same covariance matrix. 
Following this approach, it is easy to derive the following lower bound on ^(hi), 

®iAhi) = v log ( - ^ + lj (34) 

7 We note that a^o and do not depend on i, however, depends on E,=y \hj,i\ 2 - 
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where the index "g" stands for Gaussian. There are two facts that are worth mentioning about ffli ;g (hi). 
First, it is seen that lim 7 _, 0O ^'j )g (A, i ) < oo. Another point is that &i jg (hi) is an increasing function of v. 
However, setting v = u in the expression of £%i, g (hi) yields the expression of ^,fbs(^?)- Therefore, for 



all realizations of the channel gains and all ranges of 7, 



(35) 



— * 

This indicates that using &i :g (hi) as a lower bound on the achievable rate of users in the FH scheme 



provides no proof of advantage for FH over FBS. 



V. System Design 



In this section, we aim to find the optimum operation point in the FH scenario in terms of e-outage 
capacity per user. This requires finding the optimum values of {vi]f =l . For simplicity of analysis and 
fairness, we consider the fair system in which v j = v for any 1 < i < N. Therefore, the problem is 
reduced to finding the optimum value of v. As mentioned earlier in the system model, we assume the 
transmitters are not aware of the number of active users as well as the channel gains. Generally, the number 
of active users in the system is a random variable N with the probability mass function q n = Pr{A = n} 
for n > 0. We usually assume go = unless otherwise stated. 

Assuming the transmission rate of the i th user is equal to R, the outage event for this user is 



Frequency Hopping in this case. We notice that the randomness of the number of active users is also 
considered in the outage event. Therefore, the e-outage capacitjflof the FH scenario can be expressed as 




(36) 



where the subscript "FH" denotes the underlying scenario^] for which the outage is computed, e.g., 



i? FH (e) = sup {R : Pr{O i)FH (i2)} < e} . 



(37) 



The goal of this section is to find v opt given by 



v opt = arg max R FU (e). 



(38) 



V 



In the next section, the performance of the FH scheme is compared to that of the FD scheme. We distinguish the parameters of different 
scenarios by using the appropriate subscripts. 

9 We are interested in the values of e in the range e < min{q n : n G bl}. 
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We remark that v opt depends on e, 7 and {q n } n >o- 

As mentioned in the previous section, the exact expression for -Rfh(c) cannot be derived. This is due 
to the fact that a closed expression for ^fh^) is intractable. In this part, we derive lower bounds on 

— * 

-Rfh(c) using different lower bounds on M i: Fn(hi) derived in the previous section. 

— * — * — * 

Let &iyo{hi) be a typical lower bound on Mi^n{hi) for all realizations of h^ It is obvious that 

O i:F1 i( R ) C {hi : M i:lh (hi) < R}. This yields Pr{C,, FH (i?)} < Pr {h t : & ifi ,(hi) < R}, and hence, 

{R : Pr {hi : ^, lb (^) < R} < e} C [r : Pr{0 iiFH (i2)} < e}. (39) 

Defining 

i?FH,ib(e) = sup {R : Pr : M^K) < R} < e} , (40) 

we get 

i?FH,ib(e) < i?FH(e). (41) 

Based on the preceding discussion, we can derive lower bounds on -Rfh(c), namely -R F Hib( e ) an( ^ 
-R F Hib( e ) associated with the lower bounds <^ lb (/ij) and ^ lb (/ij) respectively. Consequently, we can 
obtain estimates of i> opt by maximizing -R F Hib( e ) or -^FHib( e ) over u - ^ n me following subsections, we 
separately compute -R F Hib( e ) an( ^ -^FHib( e )- 

7- Computation of -R F Hib( e ) 

We start with the following definitions. 

Definition 1 Let nGN. For c e [0, 1] and 6 > 0, a n (.; b, c) : R + — > R + w defined by 

(a , x a E B {log(l + 6(l-c B )g)}-( W -l) C log C 
a n (0; 6, c) = , (42) 

where B is a Binomial random variable with parameters (n — 1, c). 
Definition 2 Le? n e N- 2 . For 61 < 0, 6 2 > and c G [0, 1], we de/me 



/✓ ti — 1 \ m J \ 

exp ( bl 2-^-ixM j-| JJ (Mro>ro , + - n _ lfl W ■ ■ ■ d0 n _ 1; 

-i>0,-,«n-i>0 V m=lm'=l ' 

(43) 

where for each m, {0m,m' }„i'=i consists of all possible summations of m elements in the set of dummies 

{9i}^l and p m , n (c) 4 c -(l - c)"- 1 "™ 
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For example, 

^2(61, b 2 ,c) 



jf e^(b 1 2- a ^ 6 > b2 > c Xb 2 6 + l) c -(rjd6, (44) 



f exp (^-^^((Mx + l)(b 2 9 2 + l)) c(1 c) (M 2 ,i + if - 02,i V 0ld02 ' (45) 
Je 1 ,e 2 >o V / 



and 

ip 3 {bi,b 2 ,c] 

where 2 ,i = 01 + 02 by definition. 

The following Proposition offers an expression to compute -R F Hib( e ) 

Proposition 1 

[ /(1-2?)tA - 
R( FU,lb( e ) = SU P { R '■ Ql ex P + ^2 Q",1pn{bl,n, b 2, c) > 1 - 6 \ (46) 



7 



n=2 



where bi n = — , b 2 = - and c = -. 

Proof: See appendix B. ■ 
The expression given in ( [46} is quite complicated. On one hand, the multiple integrals do not have a closed 
form. On the other hand, the maximization max K lb (e) must be computed numerically. However, 
-RpHib( e ) i s tn e best lower bound on -Rfh(c) as &fi h {hi) is the best lower bound we have found on the 
achievable rate of the i th user in the FH scenario. 

2- Computation of -R F Hib( e ) 

We start with the following definition. 

Definition 3 Let n G N\{1}. For b\ < 0, b 2 > and c%, c 2 G [0, 1], we define the function <j> n (bi, b 2 , ci, c 2 ) 



as 



1 r°° 

Mh, b 2 , Cl , c 2 ) = 7 — / n - 2 exp (^(M + l)ci 2 -^ 6 ^) - 0)d0. (47) 

(n-2)!j 

Using this class of functions, the following proposition yields -R F Hib( e )- 
Proposition 2 

1 — 2v ) v\ 00 



^FH,lb( e ) = SU P { R '■ Ql eX P ( ] + ^2 Vri<Pn{bl,n, h, C 1<n , C 2 ) > 1 - 6 J> , (48) 



n=2 
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where b ljTl = ^ '—, b 2 = c 1>n = 1 - a{v,n) and c 2 = \. 

Proof: See appendix C. ■ 
Comparing the expressions for -R F Hib( e ) and -R F Hib( e )' ^ can be observed that computation of -R F Hib( e ) 
involves only one integration, while the computation of -R F Hib( e ) involves multiple integrations that are 
not tractable for many cases. To further reduce the complexity of computation, the following Corollary, 
proved in appendix D, yields another lower bound on -Rfh(^) that involves no numerical integrations. We 
denote this lower bound by -R F Hib( e )- 

Corollary 1 Let 

4tib(e) = sup Ir : f>exp (b 1>n {(n - l)b 2 + l) 1 -^ > 1 - e l . (49) 



n=l 



Then, 



R F u(e) > R<H lh (e) > 4 2 iib(e) > 4 3 H,ib(e), (50) 



where &i >n = ^ — and b 2 = J. 

Proof: See appendix D. ■ 
Fig. 1 shows the three lower bounds max„ i? F ^ lb (e) for 1 < k < 3 on max t ,i?FH(e) in a system where, 
at most, four users become active simultaneously with (qi, q 2 , q$, 54) = (0.4,0.2,0.2,0.2), u = 8 and 
7 = 100. As can be observed from this figure, max« i? F H ib( e ) an( ^ max f 4h ib( e ) 316 P re tty close to each 
other while they are within a considerable gap to max„ i2 F H ib( e )' especially for larger values of e. It is 
notable that the maximization over v is performed separately for each e. 

Having the expression for the e-outage capacity, we can find the best operational point of the system 
in terms of v, the number of selected sub-bands. For this purpose, we consider some asymptotic cases in 
terms of e and 7 and discuss the optimum value of v in these regimes. As the expression of -R F Hib( e ) * s 
not analytically tractable, we use -R F Hib( e ) f° r our analysis^ For simulation purposes, we use the most 
complex lower bound _R F ^ lb (e), which is the best bound as well. 

l0 Note that we can also use Rp 2 ^ lb (e) since it involves only one integration, however, as Rp^ lb (e) is pretty close to Rp^ lb (e), we use 
Rp% lb (e) for simplicity of analysis. 
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max T R^( E ) 
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max y R®(E) 


— <k- 
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Fig. 1. Depictions of max„ ib( e ) f° r 1 < fc < 3 in a setup where (51, 52, (ft, qi) = (0.4, 0.2, 0.2, 0.2), u — 8 and 7 = 100. 



A. Asymptotically small e 



1-2^ )v 



In this case, one can easily show that lim e ^ Rfh ib( e ) = 0- Therefore, in (|49|), we can approximate 
by — ^2_2 anc i hence, 



1 -2? ) v 



e X p P -| « 1-M (51) 

7 / 7 



By the same token, the term on the right-hand side of (49) can be approximated as 



(«) A ^In2 



Owl 



EM 1 

n=l 



7 



7 ) , (52) 



where (a) follows from the fact that the term 



\ l—a(v,n) 

/(,n l7 )4 2 *)(^ 7 + l (53) 
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does not depend on e. Using ([52]) in d49j) yields 



FH,lb 



(e) 



67 



(54) 



One can observe that the function f(v,n, 7) is concave in terms of t> for all ?i > 1. Hence, the function 
S'C 1 ') — Y^Li Qnf(v, n, 7) is concave as well. As such, the minimum of g(v) occurs either at v = 1 or 
v = u. In partiular, for t> = 1 to be the optimum value, we must have g(l) < g(u) or equivalently, 



n— 1 



n=l 



1 



a-- 



((n- 1)7 + 1) 



1- 1-: 



< 5^g„ (1 + 



rc=l 
n-1 



(n - 1)7 



E{N} 



u 



: 7, (55) 



where we have used the fact that 2 jr ( 1,n ) = 

Example 1- Assume qi = q 2 = 0.5. The condition ( |55| ) can be written as 

1\ u 



u 



u 



in— 1 



< 



1 + 7 V M 



(56) 



For each it, there is a smallest number 7^ > such that if 7 > 7„, then (56) is satisfied. For example, 
if u = 2, then 72 = 2 (3 + 2y/3). It is easily seen through simulations that j 10 < 79 < ■ ■ ■ < 73 < 72- 
However, the sequence { r y u } c ^ = i is not a decreasing sequence. In fact, one observes that the right-hand 
side of (56) tends to the increasing function in terms of 7 as u increases. On the other hand, the 



left-hand side of the same equation, i.e., the term , tends to eu as u increases. Therefore, in 



the case that both 7 and u tend to infinity, the condition (56) reduces to eu < 



1+7' 



which implies that 



7 U w lnw + lnlnu + 1. This is an increasing function in terms of u. □ 

Fig. 2 offers the curves of lb (e) in terms of v for e = 0.01, 0.05 and 0.15. The underlying network 
is characterized with qi = q 2 = 0.5, u — 10 and 7 = 20dB. It is seen that for e = 0.01, taking v — 1 
yields the best performance. As already stated in example 1, as far as q 1 = q 2 = 0.5, for sufficiently small 
e, the equation ^ (l + 10 = ^§<r — 25.8117 yields 710 = 7.052 as the minimum SNR value such 
that choosing 7 > j w results in v = 1 as the best choice. Since 7 10 < 20dB, we expect v opt = 1, which is 
confirmed in the plot. However, as the value of e increases, we are moving away from the asymptotically 
small e region and v = 1 is no longer an optimal choice. 
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Fig. 2. Sketch of Rp^ lb (e) for different values of e in a network with (31, q-x) = (0.5, 0.5), u = 10 and 7 = 20dB. 



B. Asymptotically small 7 

In this case, one can easily show that lim 7 ^ -Rfh ib( e ) = 0- Therefore, similar to the previous case, we 
can use the approximation 



1-2? « 



7 



i?ln2 

7 



Defining r = exp ^— ^-^j, one can rewrite (49) as 



^FH,lb( e ) ~ SU P \ R '■ qnT 



> 1 - e 



n=l 



where f(v,n, 7) is given in (53). As r is a decreasing function in terms of R, we can write 

7 ln(r*) 



FH.lbv / 



In 2 



where 



n=l 



(57) 



(58) 



(59) 



(60) 
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Therefore, 



11 



7 



In 2 



In ( min t* 



(61) 



It can be shown that min^ r* occurs when f(v,n, 7) takes its minimum value over v for each n. But, as 
7 <C 1, the term ( !L ^7 + l) 1 a( ~ v ' n ^ « 1 and hence, f(v,n, 7) ~ 2 jr ( 1, ' n ), which is uniformly minimized 
for all values of n by taking v = u. This gives r* = 1 — e, and hence, 



4 3 H,ib(e) = -7log(l-e) 



(62) 



This is exactly the outage capacity of a point-to-point system without interference. Therefore, in the low 
SNR regime, interference has no destructive effect on the outage capacity justifying the optimality of 



v = u. 



Fig. ^presents the plot of i?py lb (0.1) versus v for 7 = — lOdB in a system with q\ = q 2 = 0.5 and 
u = 10. It is seen that v opt — u — 10, which is expected by our analysis. Also, we obseve that the outage 
capacity is not quite sensitive to the value of v. 



Fig. 3. Sketch of lb (0.1) in a network with (qi,q 2 ) = (0.5, 0.5), 7 = -lOdB and u = 10. 
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C. Asymptotically high 7 



Recalling the expression of -RpHibl 6 ) given in (49), we have 



r( 3 ) (A 
- n 'PH,lbV fc J 



sup 



g n exp 



n=l 



jjw) (2^-1) 
7 



> 1 



(63) 



As 7 — ► 00, the term f(v, n, 7) grows polynomially with 7 for any n > 1. Assuming there exists n max such 
that Pr{iV > n max } = 0, the outage event is determined by the term with the maximum f(v, n, 7). In fact, 
since f(v, n, 7) is an increasing function in terms of n, it follows that J2™=i In exp I '-^—^ - 



vf(v 7 )(2 fl /"-l) 



. Therefore, (63) simplifies to 



^FH,lb( e ) ~ SU P { R '■ eX P 



^(^,n max , 7 ) (2^-1 )' 
7 



> 1 



which gives 



7? (3) (f) « 



f lot 



t> log 



7 m (l- 



9"- max # 



(_2-^'(^. n max) l n f 1 « ) 
j ' V <?™max/ 
/ „\"max-l 



(a) 



v log ( 1 — k^v 1 tt ' 



In 1 



(64) 



(65) 



where 



K v 



n, 



(66) 



and (a) follows from the assumption that 7 lies in the high SNR range. It is observed that the maximization 
of the above expression with respect to v is equivalent to the maximization of the term v (l — ^)™ max 1 
with respect to v, as SNR tends to infinity. This yields 



-'opt 



U 



(67) 



Fig. shows the curves of lb (0.1) versus v for different values of SNR in a system with parameters 
n max = 2, (91,(72) — (0.5,0.5) and u = 10. In general, for a sufficiently large, however finite, value of 
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SNR, one can obtain the optimum value for v by 



^opt 



arg max < v log ( 1 — n v ^ 



v \ n max — 1 



In 1 



(68) 



For example, it is easy to verify that in a system with the above parameters at 7 = 40dB, one gets 



v opt = 4, while 
7 = 40dB. 



[yj = 5. This is in agreement with the plot of lb (0.1) given in fig. j^jfor 



— »-- 


SNR=20dB 


--Q-- 


SNR=30dB 


— 1 — 


SNR=40dB 



Fig. 4. Sketch of lb (0.1) at different SNR levels in a network with n max = 2, (qi, q2) = (0.5, 0.5) and u = 10. 



Remark - In [24], we introduced a generalized version of the FH scheme (GFH scheme) where the 
number of selected sub-bands by each user changes independently from transmission to transmission. It is 
clear that the achievable rate of this generalized scheme can be higher than the case where all users only 
hop over a fixed number of frequency sub-bands. However, unlike [24] where we could analytically derive 
the optimum hopping pattern in the generalized FH scenario for the average sum-rate in the high SNR 
regime, it is not possible to find a clean mathematical formulation (or tight bounds) for outage capacity 
in the same scenario. 

VI. Comparison with other Schemes 

In this section, we compare the performance of the proposed FH scenario with that of the FD scheme 
in terms of the e-outage capacity. Frequency Division is a well-known and simple resource allocation 
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scheme that is widely used in the current wireless systems. Based on the number of existing licensees 
requesting service in the FD scenario, denoted by n des , the spectrum is primarily divided into n dcs bands^j 
and each licensed user only occupies one band upon activation. In the case that all the licensed users 
are active all the time, i.e., N = rides, this scheme results in the most efficient usage of the bandwidth. 
This makes the FD scenario superior to other schemes proposed in the literature, especially in the high 
SNR regime. However, in a practical situation, the number of concurrently active users is much smaller 
than n dcs . This makes FD highly inefficient on the heels that a considerable portion of the sub-bands is 
unused. In the simulations at the end of this section, we make an assumption that the number of active 
users N does not exceed n max , almost surely. However, n max is strictly less-^] than n des . In addition to 
FD, we also study the e-outage capacity of the FBS scenario, which is a spcial case of FH. In fact, FD 
and FBS can be considered as two extreme spectrum management schemes where the former avoids any 
interference among the users, while the latter makes all users share the same spectrum all the time. In 
the sequel, we compute Rfd(^) and -Rfbs(c). 

A. Computation of -Rfd(^) 

In the FD scenario, the spectrum is already divided into n des non-overlaping bands each containing 

^dcs 

sub-bands. Each user that becomes active occupies one of the units. As there is no interference among 
users, the outage event for the i th user can be written as 




(69) 



As hi i is a complex Gaussian random variable with variance | per dimension, \h i;i 
random variable with parameter one. Thus, 



is an exponential 




(70) 



1 1 



Each band might consist of several sub-bands. Also, it is assumed that rides divides u. 
Usually, n max < n dcs . 
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and 



i? FD (e) = sup {fl:Pr{<V D (#)}<£} 



sup < R : exp 



:i-2- 



)u 1 



> 1-e 



it 



log 1 



^-des 



™des 7 

n des 7 M 1 - e 
it 



(71) 



5. Computation of RpBs(t) 



Using (33), the following Proposition yields -Rfbs(c) 



Proposition 3 

-Rfbs(c) = sup ^ R : exp 
Proof: See appendix E. 



(l-2»)« 
7 



n=l 



(n-l)fl 



> 1 -e 



(72) 



C. Asymptotic Comparison of -Rfh(c), -Rfd( £ )> flwd i?FBs( e ) 

ij Asymptotically small e: Noting that ln(l— e) — e, we get log ^1 



n deB7 in(i-e) ^ ^ Therefore, 



win 2 



-Rfd(c) , , v 



7e 
ln2 : 



(73) 



which is the maximum achievable e-outage capacity in the underlying network for small values of e. As 
for FBS, using (72) and noting that - <C 1, we get 



-Rfbs(c) 



7 

i?ln2 



u 



71=1 ^ 



i?ln2 > 1 -e 



7 



vn=l 
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(i + 7EtT=i^ 



oo g^n-l)\ j n2 
=1 ii y 



(l + ^7) In 2 



(74) 



where in (a) we have neglected terms proportionate to R 2 . Comparing (54) and (74) reveals that R^ lb (e) is 



larger than -Rfbs(c) m the low e regime as far as (55 ) is satisfied. Furthermore, in the case that E{iV} <C -, 



the above equation implies that the FBS scheme, and consequently the FH scenario, achieves the optimal 
performance of the FD scheme. 
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2) Asymptotically small 7: Using (fTl) and ([72]), it can be realized that for the asymptotically small 7, 



R F D(e) = J R F Bs(e) = -7log(l - e), 



(75) 



which is the same value obtained in the previous section for -Rfh(^) (by setting v = u) and is the 
maximum achievable outage capacity in the network. Therefore, in this regime, spectrum division and 
spectrum sharing both achieve the optimal performance. Furthermore, the Gaussian lower bound given in 
( [34] ) implies that the FH scheme is optimal in the low SNR regime, regardless of v. However, if we use 
the proposed lower bounds (e.g., i? F 3 ^(e)), due to the presence of the term 2^>") in the expression of 
such lower bounds, taking v < u does not yield this conclusion. 



3) Asymptotically high 7: From (71), 



u , u 
log 7 H lo£ 

^■des ^des 



™dc 



U 



ln(l 



u 



log 7, 



(76) 



as 7 — > 00. Also, (72) indicates that -Rfbs(c) saturates as 7 tends to infinity and as such, FBS is highly 
inefficient in this regime. In fact, 



lim -R F Bs(e) = sup < R : } q n 2 « >l-e>. 

7— >oo ' ' 

\ n=l J 

For example, in a system where n max = 3 and (q 1 , q 2 , (&) = (0.5, 0.3, 0.2), 

lim i?FBs(e) = sup \r : 0.5 + 0.3 x 2"? + 0.2 x 2~^ > 1 - e) = u\og 



(77) 



V49 - 80e - 3 



(78) 



for any e < 0.2. 

In case of FH, ([65]) yields, 



maxi 4 3 H.ib(e) > 

V 



U 



^rriax 1 



u 



log 7 + 



u 



r?,, 



log ( — Kr u 1 In ( 1 — 

I n max I 



(79) 



For simplicity, let us assume that n max divides u. Then, 

m a.x 1 



maxR 



(3) 

FH.lb 



u 



n u 



log 7 + 



u 



log —K v In 1 — 



(80) 
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Comparing ( [76] ) and ( [80] ), we obtain 

ma.x v -Rfh(c) 



max,4 3 i lb ( e 



lim 

7^00 



RfM 



> lim 

7^00 
^dos 



1 - 



n, 



(«) n des 
> 



en r 



(81) 



where fa) follows from the fact that ( 1 



> - for any n max > 1. The above equation implies 
that n dcs > en max is a sufficient condition such that FH outperforms FD in the high SNR regime. 

Example 2- Here, we consider a practical scenario with rides = n max licensed users sharing u = n max 
sub-bands. Each user can be active with some probability p independently of other users. It is assumed 
that n max 3> 1 and p < 1 such that A = pn m3iX is a constant. Let us assume that the i th user is active. 
We are looking for a sufficient condition to guarantee a better performance in terms of the e-outage 
capacity for this user in the FH scenario compared to FD. The number of users, other than the i th user, 
which are simultaneously active together with the i th user is a Binomial random variable with parameters 
(n max — l,p). This random variable, denoted by N, can be well approximated^] by a Poisson random 
variable with parameter (n max — l)p rs A. As n max = n^es, designing the FH scenario based on this value 
is highly inefficient. In fact, for some n* > 1 where u is divisible by n*, the outage probability in the 
FH system can be written as 

Pr{C iiFH (i?)} = Pr{O itFH (R)\N < n*}Pr{N < n*} + Pt{0^ F r(R)\N > n*}Pr{N > n*} 

< Pr{O itFn (R)\N < n*} + Pr{iV > n*} 

< Pr {Oi, FR (R)\N = n*} + Pr{N > n*}. (82) 
The last line follows from the fact that for a fixed R the outage probability is an increasing function of 



the number of active users. Choosing v = ^ in (65) and selecting the transmission rate as 

R* = — log M - K^ l ~^> In ( I 



(83) 



make Pr{O i>PH ( J R*)|iV = n*} < | where 

l3 We note that N = N - 1. 



n*(n* — 1) 



Furthermore, we 
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select n* such that 



Pr{iV > n*} < -. (84) 
2 



Noting that 



Pr{iV > n*} = Pi{N > n*} 



„— A \ r, 

E < 85 > 

' n! 



n=n* 



it is shown in appendix F that for small enough e, selecting n* = — Alne guarantees Pr{iV > n*} < 



Hence, using (83), 



(a) 

RfrU) > R* 



— log ( 

n* \ 




u 7 




1 


)" 




n* 






1 


)' 




71* 




it 


log 


7, 


eA In e 



J -- ) h, (1 



!og7 + log f-Ki In (l - J 
n* \ " V 2 



log 7 

> r^!og7, (86) 



as 7 — > oo. In (86), (a) follows by the fact that setting the transmission rate at R* , we get Pr{Cj FH (i?*)} < 



e. This can be easily seen by (82) for the particular choice of n* = —Alne. 



In the FD scenario, since each active user is only allowed to utilize one frequency band, we have 

u 

# FD ( e )~ io gT . (87) 

^des 

The above equations imply that as long as n dos > — eA In e, FH is superiour to FD. This condition can be 
alternatively written as p < — -j^-. For example, if u = rides = 100 and e = 0.01, the above condition is 
satisfied for p < 0.08. □ 

D. Numerical Results 

In this section, we consider different examples to demonstrate cases where the FH scenario outperforms 
FD, i.e., max^i? FH (e) > Rfd(c) for given parameters n des , {q n }n=Ti u i e an d 7- We have no exact 
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expression for -Rfh(c)- However, we have developed the following set of lower bounds on this quantity, 

^FH(e) > 41b(e) > 41ib(e) > 4 3 iib(e), (88) 

which provide us with sufficient conditions to observe supremacy of FH over FD. -R F Hib( e ) * s me best 
lower bound, however, its computation involves multiple integrals of up to order n max — 1 for any n max . 
Computing -R F Hib( e ) on ty involves a single integral for all values of n max , whereas computation of 
i?PHib( e ) involves no integration. In the following example, we assume n max < 4. This enables us to 
use our best lower bound -R F Hib( e ) as we are a ^^ e to manipulate the double and triple integrals. In all 
numerical results, we also include the plots of the FBS scheme for the sake of comparison. 

Example 3- Let u = n des = 20, n max = 3, (91,52,93) = (0.5,0.3,0.2) and 7 = 20dB. Fig. [5] depicts 
maxt, -R F Hib( e )' -^FBs(e) an d -Rfd(c) as a function of e. It is seen that for e > 0.05 the FH scenario offers 
a better outage capacity compared to the FD scheme. Also, for e > 0.16, the FH scheme converges to the 
FBS scenario meaning that no advantage is observed by hopping over different sub-bands. Increasing the 
SNR to 30dB, fig. [6] illustrates the complete dominance of FH for e > 0.02 over FBS and FD. □ 
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Fig. 5. Comparison of FH and FD for u = nd cs = 20, n max = 3, 7 = 20dB and (qi, q2, 93) = (0.5, 0.3, 0.2). 
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Fig. 6. Comparison of FH and FD for u = rides = 20, n max = 3, 7 = 30dB and (51, q2, 93) = (0.5, 0.3, 0.2). 



VII. Conclusion 

In this paper, we considered a decentralized wireless communication network with a fixed number u 
of frequency sub-bands to be shared among N transmitter-receiver pairs. It is assumed that the number 
of users N is a random variable with a given distribution and the channel gains are quasi- static Rayleigh 
fading. The transmitters are assumed to be oblivious to the number of active users in the network as well 
as the channel gains. Moreover, the users are unaware of each other's codebooks and hence, no multiuser 
detection is possible. We considered the randomized Frequency Hopping scheme in which each transmitter 
randomly hops over v out of u sub-bands from transmission to transmission. Developing a new upper 
bound on the differential entropy of a mixed Gaussian random vector and via entropy power inequality, 
we offered three lower bounds on the e-outage capacity for each user in the proposed scheme. Asymptotic 
analysis is presented in terms of SNR and outage threshold. In the asymptotically small e regime, we 
observed that the maximum outage capacity is obtained for either v = 1 or v = u; in the asymptotically 
small SNR regime, we demonstrated that for all values of v the system achieves the optimal performance; 
for asymptotically large SNR, it is shown that t> ODt = , where n max is the maximum number of 

^ ^max 

concurrently active users in the network. We compared the outage capacity of the underlying FH scheme 
with that of the FD scenario for various setups in terms of distributions on the number of active users, 



SNR and e and showed that FH outperforms FD in many cases. 
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Appendix A; Proof of Lemma 1 

Let us consider at x 1 complex vector mixed Gaussian distribution p§(9) with different covariance 
matrices {C[}f =1 and associated probabilities {pi}f =1 given by 

L 

P&{^ = y Z / Pi9t{0,C l ), (89) 
i=i 

where g0Q) = ^^^exp [~9 H C^ l 9\. Hence, 

I p § {9) \ogp § (9)d9 = Jl < 90) 
J 1=1 

where Ji = pi J g0,Ci)logpg(O)dO for 1 < I < L. To find a proper lower bound on each J/ in this 
expression, we proceed as follows. We know that Ci = gfl t and q\ < q\ < ■ ■ ■ < q\. We have 

Ji=Pij 90, C{) log {y^p m g0, On) j dl (91) 

On the other hand, 



vm=l 



However, the term YL=i + i ttStf = Jf S ex P ~((± " tf)^*) is alwa ^ s S reater than 



Z^ m =m p , pit ■ Hence, 



i \ / L 2t l-l 



lo g lj2p m g0C m ))>lo g (p l9 0C l )) + lo g ll + £ ^ffk + E^^T I (W) 

\m=l / V 7 V m=5+l^^ m=l Pl 9t(0,Q) J 

On the other hand, the term Y! m ti ?f = £m=i ^S ex P-((i " ^f)^) is alwa ys less than 



52m=i pf p*- N° w ' we use the following inequality 



14 



which is valid for any b > and < x < a, 



log(l + 6 + x) > (1 - -) log(l + b) + - log(l + a + b). (94) 

Ci ft 



Utilizing this in the expression on the right-hand side of (93), we get: 



i (sr ra n \ \ ^ / 1 1 srPmg0,c m )\ i ^p m # f (#,c m ) 

log > Prnft^, C'm) > 1 y. TJTT^ lo S 1 + + - V 17 - log 1 + + 

(95) 

14 One may verify this using Jensen's inequality and concavity of the log(.) function. 
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where m = TL=i + i ^£ ^ "i = E™=i ?f£- Using this in yields: 



Ji>Pij 9t{e,C l )log(p l g t {e,C l ))de+ L - ^ m = lPm j log(l + /x t ) + Pm log(l + ^ + //;)■ (96) 

The first term on the right-hand side can be calculated as 



Pi J gt{0,C l )\og( Pl g t {B,C l ))de = ( Pl log Pl ) J g t {B,C x )dQ + Pl J g t {e, Q) \ogg t {9, Q)d9 

= pilogpi-pilog((7reYdetCi) 

= pi log pi - tpi Xog^egf) (97) 

where in (a) we have used the fact that the differential entropy of a t x 1 complex Gaussian vector with 
covariance matrix C\ is log ((7re)* det Thus, 

h(6) = - J p § {6)\ogp § {9)d6 

L 

= -Y, J > 

1=1 

L L 

^pi hgpi + t^pi log (negf) 
i=i i=i 

( L - ILLpA log( l + ^ + YlmtlPm lQg(1 + ^ + ^ 



< - 



1=1 

L L 



L 

t Pi log (?re^) Pi log pi 

i=i i=i 

_ ST ( L - T! m ilP m \ bg(1 + ^ + YlmtlPm lQg(1 + ^ + ^ 



L L 

t Y Pi log (^e^) Pi log Pi 

i=l z=l 



E( (pi- ) + + + Ul + ^ ) ■ 



Vl / Vl 



(98) 



Briefly, 



L 

h(Z) <tY,Pi log(vre^) + je - & 
i=i 



(99) 
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where 



Jf? = -^Pilogpi (100) 



1=1 



and 



? " = E f - E 4 lPm ) log(1 + w) + E ^ lPm log(1 + + ^ 



Sf" is a complicated function of {qi\i =1 - To simplify it, one may notice that C S" is an increasing function 
of fii. Hence, using \i x > 0, we get a lower bound on C S" , namely £f' given by 

/ — ' I// ' 

1=2 m=l 

On the other hand, using the fact that log ^ +3 ^ is a decreasing function of x, one may obtain a lower bound 
on ^ by finding an upper bound on v\ for each /. One option is v\ < ^fe ^ • Thus, we come up 
with the following lower bound on C S' 



Appendix B; Computation of -RpHit/ 6 ) 



(102) 



By (24) 



Pr {^S(^) < i?} = Pr ^ log ^ ' 7 ,„. +l)<R} (103) 

where /3 mj 7v ) = (^) m (l — ^) m and for each m, {I m)m '}l'=i consists of all possible summations 
of m elements in the set {|/^| 2 }f =ljyi , such that l£=i {^l m ,m> + l) Pm - N ^ = Ut=i(d,U+ 

We have also substituted 5fj by %^{v.,N) = aijv(2jv-i,i; H ). Since N itself is a random variable, the 
outage probability can be written as 

oo 

Pr{^i}b&) < ^} = ?i6 + E (104) 

n=2 
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where 



6 = R\N = 1} 

= Pr( v logfl + !^b^ <i? 



1-2? « 



1 — exp 



7 



(105) 



Denoting the collection of random variables {l mm >} i< m < n -i by J? n 



Pr{^ (^) < i?|JV = n} 



Pr < v log 



E < Pr < v lot 



2 -^n) 2 a„(Z_ M ;^)|^.|2 7 

2 -jr( v ,n) 2 a„(^-i,i;?,S)|/ iM | 2 7 



+ i < i? 



+i < i? 



(a) 



E < 1 — exp 



1 — E < exp 



i j_ 2 -^«.- u ,h) n n + U) 



/ 2 Jf (v,n) 



7 



1-27 « 



m=l m'=l 



_ 2 ~ a n (Z„ 



7 



n-l Ui o I v 1 

J n n + 1 

m=l m'=l 



/ 2 JT( Vl n) ( r 1 _ 2 f 

i - — 



7 t> 
f u 



(106) 



where (a) follows from the fact that after conditioning on J* n , the only random variable is \h iti \ 2 , which 
is exponentially distributed. As a result, using the definition of ip n , 

'2*M (l-2v 



r fuM = SU P { R '■ 9iexp 



1-2? v 



7 



n=2 



7 



7 v 
V u 



> 1 -e 



(107) 



Appendix C; Computation of i^Hib 



By (26), 



lib(^) =vlog I ^ — An + 1 



(108) 
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where 



E^^-IM 2 n>i 

J N = { jM . (109) 



N — 1 

Thus, following the same lines as in appendix B, we have 



Pr{^%(h t ) <R} = + (1 10) 



n=2 

where 

1-2?') v 



and 



£i = l-exp| ^ | , (111) 



(„ = I - K { exp \- ;—2- a ^'-J (lj n + l] | }. (112) 

for all n > 2. Since 2J" n ~ x\ n -iy we nave RTnW = ( n -2)\ zn ~ 2 ex P ( — ^O- Therefore, using the definition 



of 4> n , (112) can be expressed as 



/2^)(l-2fU 
l - 4> n ^ -, l 



£ n = l-0 n , -,l-a(v,n),- . (113) 

7 f M 



Hence, -RpHibl 6 ) is given by 



^FH,ib( e ) = sup < : gi exp( ) + g n 0„ , -, 1 - a(v, n), - > 1 - e 



(114) 



Appendix D; Proof of Corollary 1 

As b x < and < 2~ a ^ e ' b2 ' C2 ^ < 1 for all 6 > 0, we have 



0n(6 1 ,6 2 ,ci,c 2 ) > ? — / n - i exp(bi(b 2 9 + l) ci -9)d9. (115) 

[n — Ay. J 

One may easily check that t^^j exp(— 0)1(0 > 0) is a PDF for some nonnegative random variable 0. 
Thus, ( |115[ ) can be written as 



0n(&i, 62, ci, c 2 ) > E {exp (61(626 + 1) C1 )} . (116) 
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However, as b\ < and < c\ < 1, the function exp (61(62© + l) Cl ) is a convex function of 0. Hence, 
applying Jensen's inequality yields 

n (6i,6 2 , Cl) c 2 ) >exp(6 1 (6 2 E{e} + l) Cl ) = exp (b^n - 1)6 2 + l) Cl ), (117) 



(1 

where we have used E{6} = n — 1. Using <| 1 1 Vj) in (|48|) and noting that 61 = ^ — _J_2^(^«0 5 jj 2 = 2^ 



7 



and Ci = 1 — a(f , n), we get the desired lower bound. 



Appendix E; Computation of -R F Bs(e) 



We first compute Pr{^ )FB sOi) <R}.By (33), 



Pr{^i,FBs(/*) < = Prolog 1 + 



«(J+lL), )<i?f ' 



(118) 



where J N is defined in (109). Therefore, 



00 

Pr{^ jFBS (^) < R} = ^?nPr{^,FBsft) <i?|iV = n} 



n=l 



giPr <! wlog I 1 + \!^Aj_ ] <: /> 



+ In Prjwlog (i + ^^t|^ 7 

n=2 ^ ^ ^ « 



Jn) 



< iO ■ 



(119) 



The first term can be computed easily as 



Pr < t> log I 1 H - I < R > — 1 — exp 



7 



(120) 



For any n > 2, one can write 

\hi,i\ 2 l 



Pr < u log 1 + 



U(l + * Jn 



< R > = E <^ Pr {u log 1 + 



\Ki\ 2 i 



«(1 + W 



i7n 



(121) 



Since is an exponential random variable with parameter one, 



Pr < u log 1 + 



< R 



J n } = 1 - exp 



(l_ 2 f)(l + 27 n ) u 



7 



(122) 
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Replacing this in ( 121| ) yields 

\Ki\ 2 i 



Pr < u log 1 + 



U(l + IJn 



< R 



(1 -2«)(1 + % J n )u 




E |exp ((1 — 2^)J n 



(123) 

where we have used the fact that E{exp(tJ n )} = wz^=i as 2J n ~ xl(n-i)- Thus, -Rfbs(<0 is given by 

i , (1 — 2 u )U \ V -> _ L:// 

^FBs(e) = sup \ R : exp ( ) ^ 9n2 



> l-e>. 



7 



n=l 



Setting n* = — Alne, we have 

Pr{iV > n*} 



E 



e A 



Appendix F 



A \ n 



n=n* 



A n * 



A 



+ 



A 2 



n*\ V n* + 1 (n* + l)(n* + 2) 



+ 



< e" 



n*\ 

A n * 



E 

j=o 
1 



71* 



! 1 _ A 



n*! 1 



-A 



n* 

AeY 



1 



V2t™* V n V 1 

-A / „ \ *»' i 



n* 



y/2im* V - In e 

where (a) follows by Stirling approximation for n*\. The term (3^)™ can be written as 

_ -n*(ln(-lne)-l) 



lne 



3 A(ln(-ln e)-l)lne 



(a) 

< e 



In c 



(124) 



(125) 



39 



where (a) is valid if A(ln(— lne) — 1) > 1, which is the case for sufficiently small e. Combining this with 
the fact that I — - —K < \ gives the desired result. 
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